Sea quark polarization and semi-inclusive DIS data. 
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We investigate the potential impact of forthcoming Jefferson Lab semi-inclusive polarized deep 
inelastic scattering proton measurements in the determination of the sea quark polarization in the 
nucleon by means of a next to leading order global QCD analysis. Specifically, we estimate the 
resulting improvement in the constraints on polarized parton densities for the different flavors, 
which is found to be significant for up and strange quarks, and the correlation between remaining 
uncertainty ranges for each of the parton species. 
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U-t. '. I. INTRODUCTION 

i> : 

(N . The way in which sea quarks are polarized inside nucleons has been a persisting question ever since the spin 
, \ structure of the proton began to be unveiled by polarized DIS experiments and even today, in spite of several 
successful experimental programs, remains to a large extent unanswered. Contrary to the common belief before the 
paradigmatic EMC experiment at the end of 1980s P, the data obtained by the collaboration suggested that sea 
. quarks and gluons in the nucleon carried non negligible polarization. However this conclusion was, and has been 
C^l ■ for many years, conditional upon rather strong assumptions on isospin symmetry extended to polarized phenomena. 
, In the subsequent years, isospin symmetry itself become seriously questioned 3] and consequently the sea quark 
' polarization turned into an even more elusive question. 

Global QCD analyses including semi-inclusive measurements of polarized lepton-nucleon deep inelastic processes 
began to change this situation more recently 0, 0] , and today these data allow to constrain the extraction of polarized 
i-S^ parton densities in QCD global fits "6*1. The effectiveness of these constraints of course relies on the precision of the 
data and this is why the forthcoming generation of semi-inclusive experiments is crucial. 
O i' In a QCD global fit the uncertainty range of the resulting parton densities can be estimated by analyzing the profile 
D , of the x^-function of the fit to data against variations in the different features of the densities. This technique has 
been widely used in extractions of unpolarized parton densities _7.j and more recently has been implemented in the 
^ \ polarized case 0| providing reliable constraints on their different features, such as the net polarization carried by each 
parton flavor. 

Of course, the values each flavor polarization, or some other parameter, may take within these constraints are not 
^ independent, but become correlated. Even in the case where there are enough independent observables to extract 
■ " ' in principle all the parton densities, the uncertainties in the measurements of those observables, together with the 
theoretical uncertainties inherent in the fitting procedure, conspire against the independence of parton densities and 
results in correlations. A strong correlation between two parton densities, consequently means that neither of them is 
actually well determined. The inclusion of new and more precise data should not only reduce the uncertainty ranges 
for each flavor but also those correlations. 

In this article we investigate the impact of the inclusion of a series of semi-inclusive deep inelastic scattering 
measurements to be performed at Jefferson Lab ^ in a next to leading order (NLO) QCD global fit to all the available 
inclusive and semi-inclusive deep inelastic scattering data. In order to do this we take into account the kinematical 
coverage, statistics and level of uncertainty expected for the measurements. We also analyze the correlation between 
the uncertainty ranges of the different sea quark polarizations. As result of this analysis we find that the forthcoming 
Jefferson Lab experiment will effectively contribute to constrain the sea quark polarization in the proton. The most 
significant improvement is found in up sea quark distributions, and also with a noticeable effect for strange quarks. 
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TABLE L Inclusive and semi-inclusive data used in the fit. 
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The improvement in both distributions will be related to the inclusion in the fit of more precise data on charged 
meson electroproduction. Another interesting effect of the inclusion of these data is the softening of the correlation 
between the allowed range for the polarization of up and down sea quarks. 

In the following section, we settle definitions and conventions for the global fitting procedure and the way in which 
we study the profile of the function, we discuss the characteristics of the forthcoming semi-inclusive experiment, 
and explain how we evaluate the impact of it in a global fit. Then, we compare the results coming from the analysis 
of the set of data available at present, against those that would come from the data set enlarged with the forthcoming 
measurements, both for the individual uncertainty ranges, for the net polarization of the different flavors, and for the 
correlations between flavors. Finally, we present our conclusions. 

II. GLOBAL QCD FITS AND NEW DATA 

In the present analysis we implement the NLO QCD global fit to existing data along the lines of what was done 
in reference Q but restricting the input fragmentation functions to those of reference Q, which were shown to give 
the best fits to combined polarized data. The NLO expressions for both inclusive and semi-inclusive spin-dependent 
asymmetries and evolution equations for the parton densities can be found in 0, ^| and ^3 respectively. 

The data sets analyzed include only points with > 1 GeV^, listed in Table HJ and totaling 137, 139, and 37 
points, from proton, deuteron, and helium targets respectively, from polarized inclusive deep inelastic scattering plus 
60, 87, and 18, from proton, deuteron, and helium targets respectively from semi-inclusive deep inelastic scattering. 

The main conclusions reached in reference |^ were that using the Lagrange multiplier approach 22] as a mean to 
explore the profile of the function against different degrees of polarization in each parton flavor, definite estimates 
for the uncertainty in the net polarization of each flavor, and in the parameters of the polarized parton distributions, 
can be obtained. The overall result is a well constrained scenario where semi-inclusive data is not only consistent with 
inclusive measurements, but improves the constraining power of the fit for all the distributions, being crucial for the 
light sea quarks. 

The best fits suggest an overall picture for the quark densities at 10 GeV^ where, within uncertainties, up quarks 
are almost 100% polarized parallel to the proton, down quarks anti-parallel in a similar proportion, and sea quarks 
have a small and flavor symmetric negative polarization. The flrst moment of the gluon distribution is found to be in 
agreement with the most recent direct measurements (23j close to 0.6, constrained to be smaller than 0.8 and larger 
than -0.05 within a conservative increase in the value within a two percent range (Ax^ = 2%). 

In order to evaluate the impact of the forthcoming semi-inclusive proton data of Jefferson Lab, we included in the 
global analysis the expected values of the tt"*", tt", K'^ and K~ semi- inclusive asymmetries on a polarized proton 
target, computed with the best set of parton densities obtained in |^ with expected experimental uncertainties, as an 
additional set of points to be fitted. The projected statistical accuracies of these asymmetries are based on a total of 
225 hours of 6 GeV polarized electron beam on a polarized NH^ target. The electron beam current is assumed to be 
80 nA with a polarization of 80 %. The standard Jefferson Lab Hall C polarized NH3 target of 3 cm thickness and 80 
% polarization is assumed. The scattered electron will be detected at 30 degree with an array of lead-glass detectors 
in conjunction with a threshold gas Cherenkov counter, covering a solid angle of 210 msr. The produced hadron will 
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be detected in coincidence using the standard Hall C High Momentum Spectrometer (HMS) at 10.8 degree and a 
central momentum of 2.7 GeV/c (z^r ~ 0.5). The HMS spectrometer has a solid angle of 6 msr and a momentum 
acceptance of ±10%. 

With this enlarged set of asymmetries to be fitted, we have redone the analysis of ^ adding a detailed study of the 
correlations, and compared the resulting constraints on polarization with those of the original fit. 



III. RESULTS 



We start with the estimates for the uncertainties in the polarization of the different quark flavors. In Figure ^ we 
show the outcome of varying the of the NLO fits with the set of data available at present, in the following referred 
to as "standard fit" , and the "improved fit" , which includes the asymmetries expected to be measured by the E04-113, 
against the first moment of the respective polarized parton densities Sq at = 10 GeV^, one at a time. This is, to 
minimize 



^{Xq,aj) = X^iaj) + >^qSq{aj) q = u,u,d,d,s,g. 



(1) 



where Xq is the Lagrange multiplier associated to the polarization of a given quark flavor g, Oj are the parameters to 
be fitted, and the definition is the most simple and commonly used in fits to polarized data, namely. 
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(T,(a,)-i?.)' 



(2) 



In Equation (2), Ei is the measured value of a given observable, Ti is the corresponding theoretical estimate computed 
with a given set of parameters for the polarized parton densities, and ai is the error associated with the measurement. 
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FIG. 1: Profile of the x'^ function against parton polarizations. 
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FIG. 2: Uncertainty bands for Ap^ and Ap asymmetries coming the "improved fit" and E04-If3 expected uncertainties. 



usually taken to be the addition of the reported statistical and systematic errors in quadrature. Notice that the addi- 
tional set of asymmetries included does not contribute to when it is computed with parton densities corresponding 
to the best fit of reference 0, the same densities used to generate the asymmetries, situation that happens at the 
minima of the curves. As the distributions change in order to increase or reduce the polarization of a given flavor, 
the obtained with one or another set begin to differ. The solid lines in Figure correspond to the analysis of the 
standard set of data, while the dashed lines includes the estimated impact of future measurements. 

As expected, the most noticeable effect is in u polarization due to the maximal sensitivity of the semi-inclusive 
asymmetries on proton targets to this distribution, as discussed in j^. For the d and s distributions, the impact of the 
future measurements is comparatively suppressed by the weakness of the electric charge factor associated with these 
flavors, however there is a noticeable improvement for the s quark uncertainty near the minimum. In the previous 
analysis this distribution was mainly constrained by positivity resulting in flat distribution around the minimum 
but here shows a nice parabolic profile. Notice that both in the analysis of reference [H and in the present one, we 
are forced to assume As = As since there is not enough data yet to discern alternatives. This assumption implies a 
strong constraint on S's even though the measured asymmetries are less sensitive to this distribution than to 6u and 
dd. On the contrary, the relation between 6u and Su and the same for Sd and Sd comes from the fit. The impact on 
the gluon distribution is not significant and mainly indirect, coming as found in reference through the constraints 
on the sea quark distributions, which are now better defined. It is worth mentioning that the impact of the kaon data 
is very mild, being mostly thepion asymmetries responsible for the changes. 

As discussed in references Q , in modern extractions of parton densities it is customary to consider alternative 
sets of parton densities within an increase between 2% and 5% in x^, as a conservative estimate for the range of 
uncertainty of the global fit. In order to estimate the corresponding uncertainty range in the computation of a given 
observable, it is customary to take it as the range of variation of the observable within the alternative sets. This is 
precisely what we show in Figure |21 where we plot the uncertainty bands of and Ap~ corresponding to Ax^ = 5% 
in the improved fits as the area between the dashed lines. For comparison we include in the plots the values for these 
pion asymmetries at the kinematics of the forthcoming Jefferson Lab experiment, computed with the set of and 
which where included as "data" in the "improved fit" , together with the expected error bars. In this way, we can see 
not only the consistency of the results but also the appropriateness of the choice of Ax^ — 5% Similar uncertainty 
bands obtained for Ax^ = 5% but in the "standard fit" were found to be twice as large in 6] . 

As it was pointed out in the introduction, an important feature to keep in mind regarding the range of variation 
of the polarization of the different flavors is that they are correlated. For example, they cannot be expected to hold 
simultaneously; in consequence, for a given allowed range in x^, two or more flavors may cannot take their respective 
maximum departures from the best fit value together. For this reason it is worth while to study such correlations 
between the different fiavors, and how these correlations change with the inclusion of additional data. This can 
be done systematically generalizing Equation (1) for more than one flavor polarization, with independent Lagrange 
multipliers, scanning the proflle of the x^ function in the range of variation of them. 
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FIG. 3: Correlations between the flavor polarizations within a 5% increase in in both the standard and the improved flts. 

In Figure|31we show the allowed range of polarization within a 5% increase in in the 5u — 5d, Su — (5s, Su — 6g, 
5d — ds, 5d — dg, and 6s — Sg planes. In order to simplify the plots we have approximated the actual contours by 
ellipses, the darker ones obtained with the "improved fit" , while lighter being those coming from the "standard fit" . 
Again the most prominent effect is the shrinkage of the 5u and Ss uncertainty range. 

In these plots the correlations between the polarization of the different flavors are represented by the angles between 
the axes of the ellipses and the coordinated axes. A positive or negative 7r/4 difference would imply a maximal positive 
or negative correlation respectively, and that both polarizations are weakly constrained. This is the case, for example 
for Su and Sd in the "standard fit", situation that is corrected in the improved version. In the remaining cases, the 
axes of the ellipses are almost parallel to those of the coordinates suggesting mild correlations between the different 
pairs of flavors. 

There is however, a remaining subtle correlation between 5u and 6d and between 6u and (5s, that the enlarged set of 
asymmetries included in the fit is still not able to remove. In the first case the residual correlation is positive, while in 
the latter it is negative. The gluon polarization also seams to have negligible correlation with that of the anti-quarks, 
with only a very slight positive tendency with (5s in the standard fit, which is removed in the improved one. 

Another interesting correlation to investigate is the one between the total polarization carried by anti-quarks and 
that of the gluons. This in practice amounts to associate a Lagrange multiplier in Equation (f ) to the sum over the 
anti-quark species S'Eq — {Su + Sd + Ss) and another to the gluon polarization. In Figure^ we see that there is no 
significant correlation between ST,-q and Sg. Notice also how significantly the asymmetries expected to be measured 
by E04-II3 will help to constrain the anti-quark polarization. 
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FIG. 4: Correlation between ST/q and Sg. 



IV. CONCLUSIONS 



We have analyzed the potential impact of forthcoming semi-inclusive polarized deep inelastic scattering proton 
measurements in the determination of sea quark polarization in the nucleon by means of a next to leading order 
global QCD analysis. We find that the inclusion of this data will effectively contribute to constrain the sea quark 
polarization in the proton. The most significant improvement is found in the up anti-quark distribution, and with an 
also noticeable effect for strange anti-quarks. For down anti-quarks, the new data will have a smaller but nonnegligible 
effect. Regarding the correlations, we found that the forthcoming data will reduce the apparent correlation found 
between Su and Sd in the standard fit of reference leading to a picture where the sea quark densities and their 
uncertainties can be determine independently. 
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